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ABSTRACT 

,j  This  paper  generalizes  am  integral  representation  formula  for  eigenfunctions 

of  Sturm- Liouville  operators,  known  as  the  Volterra  transformation  operator  in  the 

I 

theory  of  the  inverse  scattering  problem,  to  higher-order  differential  operators. 
A specific  fourth-order  initial  value  problem  is  considered; 

! ' 4 d^  d 

dx  dx 

< 

(|){o)  = 1 , i|)'(o)  = 0 , (|."{o)  = -k^  , 4."'  (o)  = 0 

The  solution  for  complex  k is  expressed  as  an  inverse-Laplace-Borel  transform. 
Jump  formulae  are  obtained  relating  the  representing  kernel  directly  to  the 
coefficients  of  L.  The  result  admits  obvious  generalization  to  operators  of 
arbitrary  order. 
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SIGNIFICANCE  AND  EXPLANATION 
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The  motions  of  systems  governed  by  linear  ordinary  differential  equations  can 
often  be  developed  in  series  of  normal  modes  or  eigenfunctions;  i.e.,  solutions  of 
initial-  or  boundary-value  problems  involving  differential  equations  of  the  form 

L(J)  = 

where  L is  a linear  ordinary  differential  operator,  and  X is  a complex  number. 

In  the  theory  of  inverse  spectral  problems  — that  is,  problems  in  which  properties 
of  a differential  operator  or  system  are  to  be  inferred  from  information  about  the 
normal  mode  expansion  — detailed  information  about  the  structure  of  normal  modes 
plays  a role.  For  instance,  the  density  of  a vibrating  string  can  be  deduced  from 
observation  of  the  motion  of  a single  point  on  the  string,  which  in  turn  may  be 
interpreted  as  information  about  the  normal  modes.  In  the  solution  of  this  problem, 
an  integral  transform  of  Volterra  type,  relating  the  normal  modes  of  an  arbitrary 
string  to  trigonometric  functions  (normal  modes  of  a homogeneous  string) , plays  a 
central  role. 

The  differential  operator  which  figures  in  the  string  problem  is  of  second-order. 
Inverse  problems  involving  higher-order  operators,  however,  have  so  far  resisted 
analysis.  Among  these  is  the  problem  of  inferring  the  density  of  a vibrating  beam 
from  the  motion  of  one  point  on  the  beam,  which  again  maybe  interpreted  as  normal 
mode  information.  This  inverse  problem  involves  a differential  operator  of  fourth 
order. 

This  paper  describes  an  integral  representation  formula  for  eigenfunctions  of 
higher  order  — in  particular,  fourth  order  — operators,  which  has  many  points  in 
common  with  the  Volterra  transform  mentioned  above  in  connection  with  the  inverse 
problem  of  the  vibrating  string.  The  Volterra  transform  is  in  fact  a special  case, 
and  the  representing  kernel  is  directly  related  to  the  coefficients  of  the  operator 
involved. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  summary 
lies  with  MRC,  and  not  with  the  author  of  this  report. 
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AN  INTEGRAL  REPRESENTATION  FOR  EIGENFUNCTIONS  OF  LINEAR  ORDINARY  DIFFERENTIAL  OPERATORS 

W.  Symes 


§1.  Introduction  and  Statement  of  Results 

The  purpose  of  this  short  note  is  to  present  a type  of  integral  representation 
formula  for  the  solutions  of  linear  differential  equations  containing  a parameter  of 
the  form 

L*  = !♦  , L = .|^q^  . 

The  main  result  is  a representation  formula  for  the  solution  of  the  fourth  order  initial 
value  problem 

/ A 

X > 0 


L (x,k)  = k^i(i(x,k) 


+ r(x) 


dx  ' 


♦ (o.k)  = 1 


(o,k)  = -k  , 


dx 


— ^ (o,k)  = 0. 

dx 


(1) 


Here  q,q'  and  r are  continuous  cuid  complex-valued  on  [0,  “>  ] . 

Representation  Theorem.  Denote  by  P(x)  the  diamond  with  vertices  at 
±x,±ix.  Then  the  solution  iKx.k),  x > 0,  of  the  initial  value  problem  (1) 
admits  the  integral  representation 


1 kz 

i(x,k)  = coskx  + ® A(x,z) 


(2) 


The  kernel  A(x,z)  is  square-integrable  on  ze3P(x),  and  continuous  for  x > 0, 
ze3P(x),  z ^ ±ix;  at  the  vertices  ±ix,  A(x,z)  undergoes  a jump  given  by 


lim[A(x,  ±ix-6)  - A(x,  ±ix  + 6)  ] = --^/^dtqCt) 
6-0  ° 


(3) 


A(x,z)  enjoys  the  symmetry  properties 


A(x,  -z)  - -A(x,z)  , A(x,  z)  = -A(x,z) 


(4) 


Id 


i ^ 


: -I 

I 1 
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The  proof  is  presented  as  follows.  Several  elementary  estimates  for  the 


solution  of  (4)  are  given  in  §J.  The  main  step  in  the  proof  is  furnished  by  a Theorem 
of  B.  Ja.  Levin  on  entir'  functions  of  exponential  type,  which  is  discussed  in  §3. 

The  proof  is  completed  ir  §4. 

The  result  descrildd  above  is  clearly  a special  case  of  a general  result,  concerning 
the  initial-"eigenfunc -ion"  problem  for  a differential  expression  of  arbitrary  order  n, 
in  which  the  integration  contour  8"P  is  replaced  by  a suitable  polygon  with  n vertices. 
The  reader  will  have  no  trouble  formulating  this  more  general  result  and  supplying 
a proof  following  the  model  presented  here.  (Also,  further  jump  relations,  other  than 
(3)  are  undoubtedly  present  and  can  be  derived  by  the  method  of  §4.) 

For  n = 2,  the  contour  collapses  to  a slit,  and  the  contour  integration  is  identical 
to  integration  of  the  jump  of  the  integrand  (across  the  slit)  over  an  interval.  The 
representation  then  becomes  the  Volterra  transformation  operator  vrtiich  figures  in  the 
theory  of  the  inverse  scattering  problem  for  Sturm- Liouville  operators,  as  outlined 
in  [2].  The  jump  formula  analogous  to  (3)  is  well  known  in  that  context. 
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§2 . Some  Elementary  Estimates. 

We  consider  fourth-order  differential  expressions  on  the  half-axis  (0,“) 
of  the  form 

a / \ 

r (5) 


d^  d / d\  + 

I?  * ^ N 


with  q e C^([0,"))  , r e C*^([0,“)).  We  denote  by  the  expression  with  q = r = 0 
i.e.,  the  fourth-derivative  operator.  We  wish  to  estimate  the  solutions  of  the  initial 
value  problems 

L4)  = k%  . ♦ (0)  = , 

y U y“l 


V = 0,1,2, 3 
M = 1,2, 3, 4 , k e <t 


(6) 


To  do  this  we  view  L as  a perturbation  of  L^.  Denote  the  solutions  of  (6)  with 
in  place  of  L by  (^^(•,k).  We  have 


Jij(x,k)  = '|•(cosh  kc  + cos  kx) 

= ^(sinh  kx  + sin  kx) 


i('°(x,k) 


2k 


1 (cosh  kx  - cos  kx) 
2 


0 1 

(x,k)  = — ;r(sinh  kx  - sin  kx) 

^ 2k 

Using  variation  of  parcuneters,  we  obtain  the  following  integral  operation  for  the 
differences  i(i^(x,k)  = ij^(x,k)  - i^^(x,k): 

i|)^(x,k)  + /QG(x,y,k)i|)^(y,k)dy 


where 


X 0 

+ /^G(x,y,k)((i  w(y,)c)dy  = o 


G(x,y,k)  = — q(y)  (sinh  k(x-y)  + sin  k(x-y)) 


+ — ^ (y) (cosh  k(x-y)  - cos  k(x-y)) 


(7) 


2 dy 


2k 


+ — - r(y)  (sinh  k(x-y)  - sin  k(x-y)) 
2k^ 
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The  method  of  successive  approximations  converges  for  the  integral  equation  ( 


it  is  of  Volterra  type.  In  the  process,  we  obtain  the  estimate  for  each  x 

o 

I*  (x,k)|  < c 0 < X < X 

U o 

where  p >=  max  {|Rek|  , | Imk  | ) and  the  constant  C generally  depends  on  x 


estimates  follow  for  the  derivatives  of  ili  . 

P 


) since 
i 0 

Similar 


§3.  Some  Facts  from  the  Theory  of  Entire  Functions. 

We  collect  here  some  notions  and  facts  from  the  theory  of  entire  functions,  which 
will  be  used  in  the  next  section.  Proofs  may  be  found  in  the  monograph  of  B.  Ja.  Levin 
[5],  Chapter  I (§§  15,  16,  19,  20).  We  also  state  a slight  strengthening  of  Levin's 
generalization  of  the  Paley-Wiener  Theorem  (Appendix  I,  §3  of  [5)). 

An  entire  function  of  f is  of  exponential  type  is  there  is  A !;  0 so  that 

|f(z)|  < e^l^l  , z e 
The  (exponential)  type  of  such  a function  is  the  nuint>er 


a = lim  sup 


log  Mj(r) 


where  we  have  used  the  common  notation 


M (r)  = sv;5>  [f (re  ) | . 

^ 0ses2it 

Otherwise  put,  for  any  e > 0,  there  is  r > 0 sufficiently  large  so  that  if  |z|  > r 


|f{z)|  i 


but  there  are  numbers  z of  arbitrarily  large  modulus  for  which 


f (z)  > e 


(a-e)  z 


It  follows  that  if  an  entire  function  f obeys  an  estimate  of  the  form 


|f(z)|  < C e ' , Izl  large 

then  f is  of  exponential  type  < k.  Denote  by  E the  class  of  entire  functions  of 
exponential  type.  The  indicator  function  h^  of  an  entire  function  f of  exponential 
type  measures  the  growth  of  f in  various  directions  in  the  complex  plane.  It  is 
defined  by 


hj(0)  = lim  sup 
r -►  * 


log  I f ( re 
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The  indicator  function  enjoys  the  following  properties: 

(i)  if  f , g E E , -hen 

h^^(9)  < max{h^(e),  h^O)} 

h,  (0)  < h-(e)  + h (9) 

fg  f g 

Recall  that  the  support  function  k(9)  of  a convex  set  G in  the  plane  is  the 
infinixn  of  the  signed  distances  from  the  origin  to  half-planes  perpendicular  to  the 
ray  arg  z = 6 , which  ha 'e  no  points  in  ccsninon  with  the  interior  of  G (see  fig.  1). 

(ii)  for  f e E , is  the  support  function  of  a convex  set  , called  the 
indicator  diagram  of  f . 

(iii)  The  indicator  diagram  of  f(z)  = e is  the  singleton  set  (x) 

(iv)  If  h..(0)  S h (0)  , 0 S 0 S 2it  , then  the  indicator  of  f is  contained 

r g 

in  the  indicator  diagram  of  g . 

(v)  The  indicator  diagram  of  a sum  f + . . . + f of  functions  f . . . . f , e E is 

1 k 1 k 

contained  in  the  smallest  convex  set  containing  the  union  U 1^ 

i=l  i 


This  follows  from  (i)  . 

(vi)  Suppose  fj^,....fj^  E E , and  denote  by  L the  smallest  convex  set  containirg 
all  points  which  are  extreme  points  of  precisely  one  of  the  sets  I 

1 k 

Then  KOI  . 

1 ^ Xl^  X2^  Xk^' 

(vii)  The  indicator  diagram  of  the  exponential  sum  f(z)  = e + e +...+  e , 

with  the  X's  distinct,  is  the  polygon  with  vertices  X^^,  X^/.-.-Xj^. 

(viii)  Suppose  that  the  exponential  type  of  f is>0.  The  the  indicator  diagrcun  of 
Xz 

e f(z)  is  obtained  from  the  indicator  diagram  of  f by  translation  through 
the  vector  X. 

The  next  collection  of  facts  centers  around  the  Borel  Transform. 

If  f E E , 


f(z)  = Z 
n»0 


nl 


then  the  Borel  transform  f of  f is  defined  by  the  series 
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f(5)  = 


a> 

n=%  ^n+1 


(ix)  The  function  of  f(C)  is  holomorphic  in  the  domain  |c|  > a , where  o is 

the  exponential  type  of  f . In  particular,  f is  holomorphic  at  infinity. 

(x)  The  Borel  transform  is  additive  and  complex-homogeneous,  and  the  Borel 

X2 

transform  of  the  function  g(z)  = e is 

(xi)  Denote  by  y the  ray  arg  z = 0 . Suppose  f e e.  Then  the  integral 
0 


\ 

E 

► • • 

I . 


i 

f 


( 

i 

m 


r 


f dz  e'^^  f(z) 

'e 

converges  absolutely  and  uniformly  in  the  domain  (C:  Re  C e > h^O)  + e) 
to  f(5).  Also,  if  C is  any  circle  of  radius  = exponential  type  of 
f , then 

^ dC  "no 

The  conjugate  diagram  of  f e E is  the  convex  hull  of  the  singularities  of  f . It 
obtains  its  name  from  the  following  remarkable  theorem  of  PolyS: 

(xii)  The  conjugate  diagram  is  a function  in  E is  the  reflection  about  the  real 
axis  of  its  indicator  diagram. 

We  will  therefore  denote  the  conjugate  diagram  of  f by  1^. 

We  come  now  to  Levin's  generalized  Paley-Wiener  Theorem  ((  5 ],  Appendix  I,  §3). 
Theorem  3.1  [Levin] . Suppose  f e E,  and  suppose  that  the  conjugate  diagram  I^ 
of  f is  contained  in  the  convex  polygon  P.  Denote  by  arg  z = 9^  , j = l,...,n, 
the  rays  parallel  to  the  normals  of  the  sides  of  P,  and  suppose  that  the  functions 

-i  9 . -k ( 9 . ) X 
f^  (x)  = f(xe  ^)e  ^ 

are  square-integrable  on  the  half-axis  0 < x < where  k(0)  is  the  support 

function  of  P.  Then  the  Borel  transform  f is  holomorphic  outside  P,  and  its 

2 

boundary  values  on  the  perimeter  3P  define  a function  in  L OP).  Moreover, 
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The  proof  is  given  in  Levin's  book  ( 5 ],  Appendix  I,  §3.  We  will  also  need  the 
following  theorem,  whose  proof  is  similar  to  that  of  Theorem  3.1. 

Theori;iri  3.2.  Let  the  hypotheses  be  the  same  as  in  the  previous  theorem, 

-ie.  -k(e.)x 

except  that  the  functions  f^ (x)  » f(xe  ^)e  ^ are  assumed  to  be 

absolutely  integrable.  Then  the  Borel  transform  f assumes  continuous  boundary 
values  on  3P. 


§4.  Proof  of  the  Representation  Theorem 

To  derive  the  integral  representation,  we  first  re-interpret  the  estimates  of  §2. 
These  estimates  read 

|i|<j(x,k)|  < clkl'^eP"'  , U3(x,k)|  s clkfV’' 

where  p = max{|Re  k|,  [ira  k|}  , and 

i)<3<x,k)  = iji3^(x,k)  - y(cosh  kx  + cos  kx) 

(x,k)  = (j)  (x,k)  - — (cosh  kx  - cos  kx) 

^ 2k 

Let  (j>(x,k)  be  the  solution  to  the  initial  value  problem 

L(t>  = k‘*(ji  , 1)1(0, h)  = 1 , ().'(o,k)  = 0,  ())"(o,k)  = -k^,  I)"'  (o,k)  « 0 

and  set 


i|;(x,k)  * (J>(x,k)  - cos  kx 

= i))j^(x,k)  - k^i))3(x,k)  . 

We  obtain  the  estimate 

lii»(x,k)  I S clkp^eP’'  (8) 


in  each  finite  x-interval.  All  of  the  above  estimates  hold  for  k sufficiently  large, 
and  similar  estimates  fold  for  x-derivatives  up  to  order  three. 

Define 

h(x,0)  = X max{|cos6|,  |sine|}  , 0 < 6 £ 2it  . 


J 0 

If  we  set  k = re  , we  can  rewrite  (8)  in  the  form 

k(x,k)  < C r"^  grh(x,9)  ^ ^ 

We  can  therefore  conclude  that  iKx.k)  is  an  entire  function  of  exponential  type  ^ x , 
and  moreover  that  the  indicator  function  of  ')<(x,*)  is  bounded  above  by  h(x,9)  . 
According  to  §3,  property  (iv) , the  indicator  diagram  of  iKx,')  is  contained  in  the 
diamond  P(x)  with  vertices  ex,  ix,  which  is  the  convex  figure  of  which  h(x,e)  is 
the  support  function.  Since  ?'(x)  = P(x)  , the  Scime  is  true  of  the  conjugate  diagram 
of  ij»(x,  •)  . Moreover,  estimate  (9)  shows  that  the  function 
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is  square-integrabie  as  a fur.ction  of  r , 0 s r i “ for  0 i s 2u  . It  follows 
from  Theorem  3.1  that  the  3orel  transform  A(x,  •)  of  i|)(x,-)  is  homomorphic  in  the 
experior  of  P(x)  , assumes  square-integrabie  boundary  values  on  3P{x)  , and 

'('(x.k)  = / dz  e*^^A{x,z) 

9P(x) 

so  that  the  function  (tnx,k)  admits  the  integral  representation 
i|i(x,k)  = cos  kx  + / dz  e*^^A(x,z)  . 

2X1 

3P(x) 

We  will  now  show  that  the  kernel  A(x,2)  is  a continuous  function  of  z e 3P  (x)  , 
except  at  the  vertices  i ix  . To  do  this,  derive  frcrai  the  integral  equation  (7)  for 
= 1,2, 3, 4 , an  integral  equation  for  i(i  : 

X X 

i|)(x,k)  +/  dy  G(x,y,k)i|)(y,k)  +/  dy  G(x,y,k)cos  ky  = 0 
0 0 

where  G(x,y,k)  is  defined  in  §2,  and  satisfied  the  estimate 

|G(x,y,k)l  ^ M(y)  r ^ 

M(y)  =max{|q(y)|,  lq'(y)|,  |r(y)|) 

ie 

and  k = re 

Since  |cos  ky|  < and  h(x-y,0)  + h(y,e)  = h(x,e)  , we  obtain  for 

successive  terms  in  the  Neumann  series 


U (x,k)  = -/  dy  G(x,y,k)cos  ky 
^ 0 


0 ,(x,k)  = -/  dy  G(x,y,k)U  (y,k) 

n+1  p n 


the  estimate 


lu^(x,k)l  <(-)  e 


M-v  rh(x,0) 


M = / dyM(y) 
0 
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Write 


-X 

Jj^(x,k)  ~ ^ <3y  <3(y)(sinh  k(x-y)  + sin  k(x-y)>cos  ky  . 

Then  the  above  estimates  imply  that 

iKx.k)  = J^(x,k)  + J2(x,k) 

where  J^(x,-),  J^Cx,*)  are  entire  functions  of  exponential  type  with  indicator  diagrams 
contained  in  P(x)  , and 


J2(x,k)  « Odkl'^e"^^^*'®’) 

■ ie 

as  k = re 

It  follows  from  Theorem  3.2  that 

A(x,z)  • Jj^(x,z)  + J^(x.z) 

where  J2(x,z)  is  continuous  in  z e 3P(x).  We  therefore  need  only  show  that  Jj(x,z) 
is  continuous  except  at  z = ± ix  . To  do  this  write 

1 * 

r-  (sinh  kx  + sin  kx)  = / dt(cosh  kt  + cos  kt) 

0 


so  that 


^ X x-y 

J (x,k)  = -— / dy/  dt  q{y)(cosh  kt  + cos  kt)cos  ky 
0 0 

= -i/’'  dt/’^  dy  q(y)  + ^(-t+iy)k 

^0  0 

+ g<-t-iy)k  ^ ^i(t+y)k  ^ ^i(t-y)k  ^ ^i(-t+y)k 


Since  the  indicator  of  Jj^(x,-)  is  bounded  by  h(x,  0),  the  integral 


f dk  e J (x,k) 

0 ^ 

converges  for  Re  z sufficiently  large  to  Jj^(x,-),  the  Borel  transform  of  .J^(x,-).  It 
is  easy  to  justify  changing  the  order  of  integration.  Using  property  (x)  of  §3  we  obtain 

i 
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1 


x-t 


' I 


J (X,  z)  = - — / dt  / 

^ '*0  0 


dy  q(y)  [ 


1 


z- (t+iy) 


z- (t-iy) 


z- (-t+iy)  z+(t+iy)  z-i(t+y)  z-i (t-y) 
1 +1 


z-i (-t+y) 


z+i(t+y) 


] , 


valid  for  Re  z sufficiently  large.  An  elementary  normal  families  argument  allows  us 
to  continue  the  right  hand  side  in  the  exterior  of  P , so  that  the  eibove  representation 
is  valid  in  (t\P  . We  write  the  right  hand  side  as  a sum  of  eight  integrals.  The 
sum  of  the  first  four  defines  a function  holomorphic  in  INP  and  continuous  on 
{(t\P)u3P  , in  view  of  the  following. 

Lemma.  Suppose  that  Q(z)  is  bounded,  absolutely  integrable,  and  has  compact 
support  contained  in  the  closed  upper  half-plane.  Then  the  integral 


//  dzdz  ^ 
z— c 

Im  z>0 


converges  absolutely  and  defines  a continuous  function  of  C in  the  closed  lower 
half-plane  (C:  Im  C £ 0}  , holomorphic  for  Im  C * 0. 

The  sixth  and  seventh  integrals  can  be  shown  to  define  continuous  functions  up  to  the 
boundary  3P  : they  are  in  fact  holomorphic  in  the  complement  of  the  segment  (-ix,  ix) 
on  the  imaginary  axis,  and  have  Integrable  singularities  at  the  vertices  + ix.  The 
fifth  and  eighth  integrals  are  also  holomorphic  in  the  ccxnplement  of  [-ix,  ix]  , 
but  have  non-integrable  singularities  at  ± ix.  We  have  thus  proven  the  assertion 
regarding  the  continuity  of  A(x,z)  in  z . 

We  can  also  compute  the  jumps  of  A(x,z)  at  the  vertices  ± ix  , which  are 
precisely  the  jumps  of  the  fifth  and  eighth  integrals  ed>ove.  Ttie  computation  is 
standard;  see  for  instance  [ 1 1,  pp.  73-75.  We  obtain 
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lim  (A(x,  + ix-S)  - A(x,  + ix  +5)) 

6-»0  X 

= - j / dt  q(t) 

Since  (ji(x,k)  is  an  even  function  of  k , A(x,z)  is  an  odd 
this  follows  from  the  inversion  formula  for  the  Borel  transform 
Also,  the  reality  of  (ti(x,k)  for  k real  implies  that  A(x,z) 
z e 3P{x) . 

This  concludes  the  proof  of  the  representation  theorem. 


function  of  z : 

§ 6,  property  (ix) ) 
= - A(x,z)  for 
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